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Abstract 

The two-dimensional S = 1/2 asymmetric Heisenberg Mattis model is investigated with the 
exact diagonalization of finite clusters. The Neel order parameter and the spin glass order 
parameter can be smoothly extrapolated to the thermodynamic limit in the antiferromagnetic 
region, as in the pure Heisenberg antiferromagnet. The critical concentration of the Neel phase 
is consistent with that of the two-dimensional Ising Mattis model, and the spin glass order 
parameter increases monotonously as the ferro-bond concentration increases. These facts sug- 
gest that quantum fluctuation does not play an essential role in two-dimensional non-frustrated 
random spin systems. 

KEYWORDS: quantum spin system, ground state, randomness, Mattis model, Neel order, 
spin glass order 

1 Introduction 

In recent studies on spin systems, effects of quantum fluctuation and randomness have been 
investigated from various aspects. In considering recent development of computer facilities, it is 
high time to study the coexistence of these two nontrivial effects. In fact, numerous studies have 
already been made for the inflnite-range Ising model in a transverse fleld[0 and the inflnite- 
range quantum XY model. Quite recently, analytic and numerical studies of the short-range 
random Heisenberg model began to be attempted in one-^, ^, |^ and two-dimensional [|^, 
systems. However, all these works are motivated by the question how ground-state properties of 
quantum spin systems are affected by random perturbations, and the behavior of characteristic 
quantities in random spin systems such as the spin glass (SG) order has not been considered 
until present. 
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As a first step to quantum SG problems, we investigate the two-dimensional S = 1/2 
asymmetric Heisenberg Mattis model on a square lattice described by the following Hamiltonian, 

H = JijSi ■ Sj , (1.1) 

where (ij) denotes all the nearest-neighbor pairs, with 

PiJ,f,p) = p5{J,, -J) + (1 -p)5(J,, + J) , (1.2) 

Jij = JTiTj , Ti = ±1 . (1.3) 

Since frustration does not exist ||^ in the present model, the coexistence effect of quantum 
fluctuation and randomness is expected to be displayed clearly. As is well-known, the free 
energy of the following Ising Mattis model, 

n = -j2j,^s^s] , s = i, (1.4) 

(ij) 

with the conditions (|1.2|) and (|1.3|) , is equivalent |@] to that of the corresponding pure Ising 
model through the following gauge transformation of the spin variable {S^}, 

St ^ Sin . (1.5) 

On the other hand, such equivalence does not hold in the quantum Heisenberg Mattis model, 
and ground-state properties of this model are nontrivial. In the present paper, we investigate 
the behavior of the Neel order parameter and the SG order parameter. 

In §2, properties of the asymmetric Ising Mattis model are briefly reviewed, following Ozeki's 
paper. IP His argument is extended to the asymmetric quantum XY Mattis model. In §3, 
the finite-size correction of the SG order parameters is investigated by means of the exact 
diagonalization of the two-dimensional = 1/2 XY model. In §4, ground-state properties of 
the two-dimensional S" = 1/2 Heisenberg Mattis model are investigated similarly. The critical 
concentration of the Neel phase and the p-dependence of the SG order parameter are 
calculated. In §5, these descriptions are summarized. 



2 Review of the Asymmetric Mattis Model 

Even in the classical Ising-spin case, the Mattis model p| shows nontrivial phase transitions [§] 
when asymmetric bond distribution {p ^ 0.5) is considered. Namely, the phase diagram of the 
Ising Mattis model ( |1.4| ) is shown in Fig. 0, and this phase diagram can be obtained analyt- 
ically from that of the corresponding pure Ising model. This mechanism can be understood 
through the behavior of correlation functions. Following Ozeki's paper, let us consider the 
ferromagnetic correlation function, 

g{r-p,T)^[{S^,S:)l , (2-6) 
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and the SG correlation function, 



^(r;p,T)^K5o^5;)^]r. (2.7) 

Here (■ ■ ■) denotes the thermal average, and [■ ■ -jj, represents the random average based on the 
conditions (|1.2| ) and ( p..3|) . After the gauge transformation (|1.5|) , these two correlation functions 
can be expressed as 

g{r-p,T) = (5o^5;)p,re X [rox^ , (2.8) 
~g{r;p,T) = {S^S:)^,, x [{ToTr)\ 

= {SqS^)p^j.^ , (2.9) 

where (■ ■ ■)pure denotes the thermal average in the pure Ising model {p = I case of the model 
( |1.4| )). Then, the structure of the phase diagram given in Fig. ^ can be explained clearly. 
The para-ferro (or the para-Mattis SG) phase transition originates from the singularity of the 
correlation function (S'oS'^)purc, and the ferro-Mattis SG phase transition originates from that 
of [toTj-Jj,. Since the latter correlation function is independent of temperature, the ferro-Mattis 
SG phase transition occurs even in the ground state. Moreover, this random average can be 
related^ with the thermal average (tqT^) defined in the following Ising model, 

^fW = -E^^^J for p> 1/2, (2.10) 

HafW = for P< 1/2. (2.11) 

In this argument, the difference between the micro-canonical ensemble (the original bond dis- 
tribution based on the conditions ( |1.2|) and (|1.3| )) and the canonical ensemble (with respect 
to the Ising models ( 2.101 ) and ( 2.11|) ). Then, the critical concentration p^^^^ of this model is 



related]^ with the critical temperature Tc of the Ising models ( p.lO| ) and (|2.11|) , as 



eF(r,) = l-2pf, (2.12) 
eAF(Tc) = 2pf-l, (2.13) 

where ep and caf represent the energies per bond of Hf{t} and T^afI''"}, respectively. There- 
fore, these two phase transitions at T = Tc or p = p^^^^ belong to the Ising universality class. 
Especially, in the two-dimensional Mattis model on a square lattice, the critical temperature 
and the critical concentration are given0 by Tc = 2.269 ■ ■ ■ and = 1 — p^^ = 0.853 ■ ■ ■, 
respectively. 

Ozeki's argument is valid even in quantum spin systems, as long as they can be transformed 
to pure systems through similar gauge transformations. ||10|| Namely, the quantum XY Mattis 
model described by the following Hamiltonian, 

n = -Y.Ms:sj + sfs]) , (2.14) 



with the conditions ( p..2|) and (|1.3| ), can be treated similarly through the following transforma- 
tions, 







sf - 


SfTi 



(2.15) 
(2.16) 

S- - St. (2.17) 

Since the existence of the long-range order has already been proved in the pure two-dimensional 
S = 1/2 XY model at the ground state, |]TT[ the ground-state phase diagram of the model 



( |2.14| ) should be equivalent to that of the two-dimensional Ising Mattis model: The critical 
concentrations are given by = 1 —p^^ = 0.853 ■ ■ ■, and the value of the SG order parameter 
is independent of p. 



3 Finite-Size Correction of the Spin-Glass Order Param- 
eter 

The spin-wave theory tells that the asymptotic form of the two-point correlation function of 
the two-dimensional antiferromagnetic Heisenberg model {p = case of the model ( |1 . 1| ) ) at the 
ground state is given |T2|, |1^] by 

0,. S,),r^ ( mo + ^J—] , (3.18) 

with Tij = \ri — rj\, and (■ ■ •)g represents the ground-state average. Thus, the Neel order 
parameter mst{N) defined in 



^ ■ S,)l , (3.19) 



and the SG order parameter msg{N) defined in 



<iN)^^^J:[{S^■S,fl , (3.20) 



are expected to be scaled as 



ml{N) ~ + const, x A^-^/^ ^ o{N-^) , (3.21) 
mUN) ~ + const, x A^-^/^ + 0{N-^) . (3.22) 



These formulas mean mst{oo) = mo and msg(oo) = mp, and the scaling form ( |3.21 



was 



confirmed by quantum Monte Carlo simulations p!4|, 15| of this model. Although these scaling 



forms are not trivial when randomness is introduced, essential behavior is expected to be similar 



in the Neel phase, as in the diluted quantum antiferromagnet. ITq, \l% On the other hand, the 
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Neel order does not exist in the Mattis SG phase {p ~ 0.5), and the asymptotic form of the 
two-point correlation function is supposed to be 

(5. -5,)^- (-1)1-^1 exp(^-i-^ . (3.23) 

Then, the scahng forms of the order parameters are given by 

ml{N) ~ const. X A^-i + 0(iV"2) , (3.24) 
ml^{N) ~ const, x N'^ + 0{N-^) . (3.25) 

Namely, the existence of the SG order parameter in the Mattis SG phase cannot be derived 
from the asymptotic form of the two-point correlation function. Therefore, the scaling form 
( 3.22 ) obtained from ( p.l8 ) is also questionable. In fact, when m'^ (N) of the antiferromagnetic 



Heisenberg model is plotted versus N'^f^ as predicted in ( p.22| ) (Fig. ||), this order parameter 
does not seem to exist in the thermodynamic limit, which is not consistent with the existence 
of the Neel order parameter. 

Then, the size dependence of the SG order parameter can only be argued numerically at 
present. For this purpose, the two-dimensional S = 1/2 XY model is considered. According 



to spin- wave analyses and quantum Monte Carlo simulations of this model, fl^ the size de- 
pendence of the Neel order parameter is equivalent to that of the two-dimensional S = 1/2 
antiferromagnetic Heisenberg model. This fact is due to the similarity of the structure of low- 
energy excitations in these two models, and the size dependence of the SG order parameter of 
the XY model is also expected to be as such. Moreover, even if the Mattis-type randomness is 
introduced, the value of the SG order parameter does not vary in this model, as explained in 
the previous section. Since the x, y-ax.es and the 2;-axis are not equivalent in the XY model, we 
should consider two different SG order parameters defined in 



iV2 



= (3-26) 

and 



«Wr-42E[(^r^|)l ■ (3-27) 
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We diagonalize Oitmaa-Betts-type clusters [0 up to the 20-spin one, and calculate these order 
parameters. The quantity m^g should vanish in the thermodynamic limit, and m^^ should 
remain finite. In addition, the size dependence of these two order parameters is expected to 
be similar. All these conditions are satisfied when these order parameters are plotted versus 
A^~^, as shown in Figs. ^ and |^. It is true that we cannot rule out the possibility of ~ A^~^ 
(6' ~ 1, but not 1) dependence only from these numerical data, but such a fractional exponent 
can only appear when the system is critical. Thus, we assume the scaling form of the SG order 
parameter as 

ml^{N) ~ m2g(oo) + const, x A^^^ + 0{N-'^) . (3.28) 
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4 Numerical Results of the Heisenberg Mattis Model 



As in the previous section, we diagonalize Oitmaa-Betts-type clusters up to the 20-spin one. 
Note that this cluster size is not so "small" in the sense of numerical calculation, because spatial 
symmetries cannot be used in the Hamiltonian-matrix diagonalization in random spin systems. 
In compensation for such limitation of cluster sizes, we average all the bond configurations. 
Since spatial symmetries can be applied to this summation, such calculations are within the 
reach of recent super computers, though still tedious. Then, the present results are "exact", 
not only in the sense of Hamiltonian diagonalization, but also in the sense of random averaging 
in finite clusters. 

In an A^-spin Oitmaa-Betts-type cluster, the number of bonds are 2N , and the number 
of ferromagnetic bonds which satisfy the conditions ( |1.2| ) and ( |1.3D are limited to = 
0,4,6, ■■■ ,2N — 4,2N. Then, in order to calculate a physical quantity Q for an arbitrary 
concentration p, we should extrapolate the values of Q at p = Np/2N. In the present pa- 
per, we sum up all the bond configurations on the basis of the canonical distribution^ of the 
corresponding A^-spin Ising model. Namely, using the "temperature" T{p) determined by 

eF{T{p)) = l-2p for j9 > 1/2 , (4.29) 

eAF{T{p)) = 2p-l for p< 1/2, (4.30) 

we have 

,_ Ep^^(p)g(p)e-^^(^)/^BT(p) 

Epw{p)e-^^(py^Bnp) ' ^^"^^^ 

where w{p) denotes the number of bond configurations in which the ferro-bond concentration 
is equal to p, Q{p) represents the partial random average restricted to the samples in which the 
ferro-bond concentration is equal to p, and Eg{p) is given by 

Eg{p) = 2N{l-2p) for p> 1/2, (4.32) 
Eg{p) = 2N{2p-l) for p< 1/2. (4.33) 

At first, we estimate the Neel-Mattis SG phase boundary by means of the least-squares 
fitting of the Neel order parameter m^^. In Fig. ^ m'^^^ for the 8, 10, 16, 18, 20-spin clusters are 
plotted versus 

jY-i/2 p ^ Q Q Q i^^ Q 2. This fi gure shows that m^^ vanishes at p ~ 0.15. 
Then, we vary the value of p more precisely, and obtain the following estimate of the critical 
concentration of this phase boundary, 

= 0.140 ±0.016 . (4.34) 

This error bar means that the estimate of rrist does not coincide with zero for p < 0.124 
and p > 0.156, within the standard deviation in the fitting. This estimate is consistent with 
that of the Ising Mattis model, p^^ ~ 0.147. Although clusters are still small, the fitting 
curves are linear enough, as in the pure quantum antiferromagnet. The reason of such a good 
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scaling behavior may be that wavefunctions are not locahzed in antiferromagnetic random spin 
systems owing to quantum fluctuation. On the other hand, similar analysis is not possible for 
the ferromagnetic order parameter, because quantum fluctuation is small in the ferromagnetic 
region. 

For p > p^^, the extrapolated value of m'^^ is negative, which simply means that the scaling 
form (|3.21|) based on the existence of the Neel order is not correct. The expected scaling form in 
this region is ( |3.24| ), and this form is satisfied in the region p ^ 0.5. The linearity of the fitting 



curve just above the critical concentration is expected to be a finite-size effect. Moreover, the 
size dependence of the Neel order parameter can be fractional at p = p^^ as m'^^{N) ~ A^~^, 
with 1/2 < 6 < 1. These facts show that the estimate ( [4.34|) may be underestimate. However, 



as long as p < p^^, the size dependence ( p.21|) holds no matter how small mgt is, and the 
present naive fitting can be justified practically. In fact, the critical concentration of the 
diluted antiferromagnetic Heisenberg model was analyzed similarly, [l^ Furthermore, it is quite 



improbable that the critical concentration is increased by quantum fluctuation. Therefore, the 
good coincidence of the estimate ( [4.34| ) and the classical value p^^ ^ 0.147 can be regarded as 
the evidence of the smallness of such finite-size effects. 

Next, we investigate the stability of the Mattis SG phase against the coexistence of random- 
ness and quantum fluctuation. In Fig. ^ the SG order parameter m^g for the 8, 10, 16, 18, 20-spin 
clusters are plotted versus 1/A^ for p = 0.0, 0.3, 0.5, 0.6, 0.7. The extrapolated value of this order 
parameter increases monotonously, as p increases. This fact shows that this order parameter is 
reduced only by quantum fluctuation, in spite of the coexistence of randomness. Although the 
assumed size dependence ( p.28| ) is not based on a strict theoretical background, such monotony 
of the estimate is independent of scaling forms. Therefore, as long as the Neel order exists at 
p = 0, the above result will not be changed. 



5 Summary and Discussions 

The present numerical study suggests that the ground-state phase diagram of the two-dimensional 
S* = 1/2 asymmetric Heisenberg Mattis model is equivalent to that of the two-dimensional Ising 
Mattis model. Namely, the values of the antiferromagnetic critical concentration p^^ are consis- 
tent within the fitting error in these two models, and the SG order parameter in the Heisenberg 
Mattis model increases monotonously as p increases. These facts show that magnetic order- 
ings in the two-dimensional Heisenberg Mattis model is mainly controlled by randomness, and 
quantum fluctuation only reduces the size of spins. If we can identify quantum fluctuation with 
thermal fluctuation, this result can be understood clearly, because the Neel-Mattis SG phase 
boundary is independent of temperature (i.e. vertical) in the Ising Mattis model. If this picture 
can be justified, similar quantum-classical relationship is expected in the two-dimensional ±J 



SG model, because the boundary of the Neel phase has been proved vertical [|19|, |2^ in the 
classical case. 
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The present study also shows that ground-state properties of random quantum spin systems 
can be evaluated from relatively small clusters in the antiferromagnetic region. This fact is 
quite encouraging for further studies on ground-state properties of quantum spin glasses. Some 
studies in this direction |21|, |23] are now in progress. 
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Figure Captions 



Fig. 1: Phase diagram of the Ising Mattis model. denotes the critical point of the corre- 
sponding pure Ising model, and is related with T^.. [(5'^)]r = and [(>S'^)^]i. = in the para 
phase, [{S^)]r ^ and [{S^f], in the ferromagnetic phase, and [{S^)]r = and [{S^)\ ^ 
in the Mattis spin glass phase. 

Fig. 2: The spin-glass order parameter of the two-dimensional S = 1/2 antifcrromagnetic 
Heisenberg model on a square lattice is plotted versus A^^^/^. The number represents the size 
of clusters, and the straight line is drawn by the least-squares fitting. 

Fig. 3: The spin-glass order parameter (^-component) of the two-dimensional S — 1/2 XY 
model on a square lattice is plotted versus 1/N. The number represents the size of clusters, 
and the straight line is drawn by the least-squares fitting. 

Fig. 4: The spin-glass order parameter (x-component) of the two-dimensional S = 1/2 XY 
model on a square lattice is plotted versus The number represents the size of clusters, 

and the straight line is drawn by the least-squares fitting. 

Fig. 5: The Neel order parameter of the two-dimensional S = 1/2 Heisenberg Mattis model is 
plotted versus N^^/'^ for p = 0.0, 0.1, 0.15, 0.2. The straight lines are drawn by the least-squares 
fitting. 

Fig. 6: The spin-glass order parameter of the two-dimensional S — 1/2 Heisenberg Mattis 
model is plotted versus 1/N for p — 0.0, 0.3, 0.5, 0.6, 0.7. The straight fines are drawn by the 
least-squares fitting. 
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Figure 4: 
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